Modified Harris Hawks optimization for the 3E feasibility assessment of a hybrid renewable energy system

The off-grid Hybrid Renewable Energy Systems (HRES) demonstrate great potential to be sustainable and economically feasible options to meet the growing energy needs and counter the depletion of conventional energy sources. Therefore, it is crucial to optimize the size of HRES components to assess system cost and dependability. This paper presents the optimal sizing of HRES to provide a very cost-effective and efficient solution for supplying power to a rural region. This study develops a PV-Wind-Battery-DG system with an objective of 3E analysis which includes Energy, Economic, and Environmental CO2 emissions. Indispensable parameters like technical parameters (Loss of Power Supply Probability, Renewable factor, PV fraction, and Wind fraction) and social factor (Human Developing Index) are evaluated to show the proposed modified Harris Hawks Optimization (mHHO) algorithm’s merits over the existing algorithms. To achieve the objectives, the proposed mHHO algorithm uses nine distinct operators to obtain simultaneous optimization. Furthermore, the performance of mHHO is evaluated by using the CEC 2019 test suite and the most optimal mHHO is chosen for sizing and 3E analysis of HRES. The findings demonstrate that the mHHO has achieved optimized values for Cost of Energy (COE), Net Present Cost (NPC), and Annualized System Cost (ASC) with the lowest values being 0.14130 $/kWh, 1,649,900$, and 1,16,090$/year respectively. The reduction in COE value using the proposed mHHO approach is 0.49% in comparison with most of the other MH-algorithms. Additionally, the system primarily relies on renewable sources, with diesel usage accounting for only 0.03% of power generation. Overall, this study effectively addresses the challenge of performing a 3E analysis with mHHO algorithm which exhibits excellent convergence and is capable of producing high-quality outcomes in the design of HRES. The mHHO algorithm attains optimal economic efficiency while simultaneously minimizing the impact on the environment and maintaining a high human development index.

At present there is a growing global inclination towards the advancement of Renewable Energy Sources (RESs).Consequently, there has been an expanding imperative to cultivate sustainable and eco-friendly energy sources commonly referred to as RESs.Solar, wind, small hydro, cogeneration bagasse, and biomass, are just a few of the sources of renewable energy that have great potential to be used to power developing economies.For example on March 31, 2021, it was predicted that India has a total renewable power generating capacity of 14,90,727 MW 3 and the full overview of the contributions provided by different RES in India is shown in Fig. 1.
Each specific RES exhibits a unique set of advantages and disadvantages.Solar and wind energy as sources of renewable energy have become known for their cleanliness, affordability, and environmental sustainability.It is important to note that their effectiveness is contingent upon various factors such as weather conditions, climatic patterns, and specific geographical regions.Despite the intermittent and unpredictable nature of solar and wind sources, which can be attributed to abrupt changes in solar radiation as well as the speed of the wind, they can serve as viable substitutes for fossil fuels in power generation 4 .Due to their unpredictable nature, RESs that rely on a single source, like solar and wind power, are unable to provide the year-round need of electricity.However, the hybridization of solar Photovoltaic (PV) and Wind Turbine (WT) technologies have the potential to fulfill demand for electricity under the condition of adequate solar irradiation and wind speeds.As said already, the solar and wind sources are intermittent, thus there exists a requirement for an Energy Storage System (ESS).Numerous technologies for storing energy have been implemented in various applications, encompassing batteries, hydrogen, compressed air, flywheels, pumped hydro storage, and gravity energy.Batteries are extensively utilized as a short-term storage medium owing to their affordability, ease of upkeep or maintenance, and consistent operational capabilities.Green hydrogen possesses several notable benefits including its ability to store energy over extended periods, and its capacity to mitigate local pollutants in an environmentally friendly manner.Also, the specialty of Green Hydrogen is the complete absence of CO 2 emissions.However, it is important to acknowledge that green hydrogen is accompanied by a comparatively higher cost.At present, the water electrolysis process stands as the most ecologically sustainable method for hydrogen production with a drawback of substantial electricity consumption 5 .The utilization of hydrogen fuel cells facilitates the conversion of the energy stored within hydrogen into electrical energy.
It has become increasingly clear that a Hybrid Energy System (HES) that incorporates mixed generation and suitable storage is necessary to meet the requirements for a consistent, sustainable, and dependable power supply.As a result, determining the optimal sizing for these systems has become a highly desired research topic worldwide.By implementing appropriate component sizing techniques and adopting an effective power management strategy or energy dispatch strategy, it is possible to establish a Hybrid Renewable Energy System (HRES) that offers a sustainable, economically viable, and dependable solution for fulfilling the energy requirements of developing nations.A key consideration in the design of HRES is the sizing of its components.Solar and Wind energy sources exhibit discontinuity, dependence on specific locations, and significant variability 6 .Consequently, the economic viability of stand-alone energy systems is compromised due to the inability to properly match power generation with energy demand fluctuations 7 .HRES which utilize solar PV, wind technologies, ESSs, and Diesel Generators (DGs) have the potential to serve as viable energy sources 8 .
The optimization and efficient operation of HRES will contribute to the enhanced dependability, reliability, and control, along with the efficiency of the power network.Efficient management of these components is crucial for enhancing the HRES operation or performance.The implementation strategy of HRES is determined by technical, economic, and environmental indicators within the network 9 .In recent times, there has been a significant focus on the identification of the most suitable HRES configuration and capacity.The HRES is established purely by taking into consideration the presence of RESs and annual power demands specific to a given location.The determination of a suitable size for the HRES through the application of techno-economic evaluation has been a prominent area of research for the past quarter-century or over 25 years.
Literature explores the optimum size of off-grid systems to reduce the Net Present Cost (NPC), Cost of Energy (COE) or Levelized Cost of Energy (LCOE), and Loss of Power Supply Probability (LPSP) while maximizing the renewable factor (RF) 10 .SK.A. Shezan et al. performed a practical evaluation of an independent WT/PV/ DG/battery HES designed for rural areas.The HOMER software tool was employed to model and enhance the efficiency of the PV/WT/DG hybrid energy system using real-time field data of solar radiation as well as wind speed pertaining to the given area.The simulations using HOMER software usually confirms that the framework is feasible in relation to NPC analysis and the idea of decreasing CO 2 emissions 11 .The various metrics including Fig. 1.Different RES contributions in India.
Vol:.( 1234567890 Investigated the design and performance of a desalination system that was powered by solar energy.Although the model identifies substantial economic and environmental advantages, its long-term efficacy cannot be confirmed due to the absence of real-world implementation data.Furthermore, it relies on HOMER software, which might limit the complexity of real-world system modeling 40 Gas Emission (GHGE) savings, Human Development Index (HDI) and maximum reliability of the system.The Harris Hawks Optimization (HHO) algorithm proposed by Ali Asghar Heidari et al. 50has attracted attention in the field of optimization challenges.According to the literature, the HHO algorithm is applied to the sizing, optimization, and design of autonomous MG, and various applications have been proposed recently 51 .Further renewable energy-based DG unit sizing and planning for the IEEE bus system are also analyzed 52 .Gauri Sahoo et al. presented an approach including the utilization of Modified HHO to develop a Fractional-Order Fuzzy Proportional-Integral-Derivative (PID) controller.The controller was aimed to regulate the frequency of a multi-MG system 53 .Mohammed Kharrich et al. developed an equilibrium optimized for the techno-economic analysis of MGs and the results are compared with the HHO algorithm 54 .Ch Srivardhan Kumar et al. provided a control approach that offers a reduced level of complexity for voltage compensation to boost power quality inside the distribution systems.The control technique proposed by them utilized a customized power device known as the Dynamic Voltage Restorer (DVR), which was implemented with various inverter topologies with the application of the HHO algorithm 55 .Zhang et al. proposed a method for phase partitioning using HHO with a hard sequential constraint.It was proposed to find the optimal phase partitioning results within a specific target phase number in the inner loop.Additionally, it automatically finds the optimal phase number in the outer loop by balancing the trade-off between modeling complexity as well as partition execution 56 .
The fundamental concept of this HHO technique relies on the hunting strategy employed by the Harries Hawk in capturing its prey.The HHO algorithm is a relatively new and highly promising optimization methodology for several reasons.The HHO theory is based on a population approach and does not rely on the computation of costly partial derivatives.As a result, it is classified as a gradient or derivative-free optimization method 57 .The aforementioned benefit enables the utilization of HHO in the optimization of any problem, provided that the problem is well formulated.It is recognized as a global optimization algorithm due to its incorporation of both exploitation as well as exploration phases.However, the optimal convergence curve as well as local minima are two challenges reported in the literature.To mitigate this challenge a modified version of the HHO (m-HHO) algorithm is proposed in this work.
The research proposed in this work utilizes a mHHO method to identify the optimal size of the HRES system.This is the first and innovative approach where HHO is modified with nine different operators using inertia weight strategies.The effectiveness of these improved algorithms has been tested by utilizing CEC 2019 test functions 58 .In addition, the statistical importance of this algorithm is determined by using a statistical test known as the Friedman rank (f-rank) test 59 .This statistical test is non-parametric in nature.The f-rank test involves assigning a distinct rank to each algorithm being compared based on its performance.The total f-rank is then determined by considering the results across all test functions.Box plots and statistical test studies have also been performed to check the accuracy as well as the robustness of the proposed algorithm.From the acquired results best mHHO is selected for HRES design and the results of the design are compared with some well-known MH-algorithms such as Differential Evolution (DE) algorithm 60 , Enhanced Whale Optimization Algorithm (EWOA) 61 , Particle Swarm Optimization (PSO) 62 , jDE100 63 , Novel Bat Algorithm (NBA) 64 and Sin Cosine Algorithm (SCA) 65 .Also, the f-rank test and Wilcoxon test (p-rank) are carried out to verify the effectiveness and robustness of the mHHO algorithm for picking the optimal mHHO to design the hybrid system.
From the literature, it is noted that the environmental, economic, and ecological (3E) assessment along with social and technological elements on optimal sizing of HES has not been reported in recent studies.Also, it is found that the use of a combination of COE-LPSP-HDI to establish the optimal size of HRES is not often seen in the current literature.Therefore, this study specifically focuses on several objectives to perform 3E evolution and a comparative study of the proposed HES system using MATLAB simulation.Here the sizing and 3E analysis of the HRES system is presented as an objective function defined as the minimization of COE, NPC, ASC, CO 2 emission, and LPSP thereby improving PVF, WTF, RF, and HDI.Also case study of the HRES system consists of PV-WTG-BESS-DG configuration.The contributions of the work are summarized as follows: • A modified version of HHO known as the Modified HHO (mHHO) algorithm is presented in this paper and this is the first such attempt in the HRES design utilizing mHHO for rural electrification.The subsequent sections of the document are organized in the following manner: "Mathematical modeling" of the paper outlines a mathematical model of the various components of an HRES."Case study" emphasizes the details of the selected location for a case study."Methodology" describes the methodology which includes a power management strategy and the formulation of an optimization problem, which further consists of two components: objective functions and system constraints."Proposed optimization technique" demonstrates the basic HHO algorithm and proposed optimization technique mHHO algorithm."RESULT AND DISCUSSION" www.nature.com/scientificreports/ is divided into two subsections where the first section analyzes and discusses the simulation outcomes of the mHHO algorithm and the second section addresses the simulation results of HRES using mHHO.Finally, "Conclusion" offers few conclusive findings or remarks.Figure 2 illustrates the detailed structure of the proposed work.

Mathematical modeling
Hybrid Power Systems (HPS) integrate a minimum of one RES (such as wind, solar, or hydropower) with one or more traditional sources of energy, such as DGs.Due to its reliance on atmospheric conditions as already pointed out, the production of electricity from an RES is unpredictable.Therefore, combining this with conventional power sources and/or battery storage would provide reliable and constant electrical power and in this work such a system is considered as HRES.A HRES is typically made up of generation, distribution, and demand subsystems, which may all change significantly based on a number of factors.Some of the factors are the accessibility of renewable resources, the intended services to be offered, and the demand subsystem.The aforementioned parameters exert a significant influence on the decision-making process, and consequently on both the cost as well as dependability of the system 66 .
The design of HRES includes several elements, including a solar PV power generation system, WT power generation system, Power Converters (PC), Battery Storage System (BSS), and DG system which are depicted in Fig. 3. Distributed generation is typically utilized as a final option and is exclusively started when other energy sources are incapable of satisfying the load demand at a specific moment.RESs are expected to offer a consistent supply of electrical power throughout some interval which depends on various factors and optimization of those factors is of utmost interest in this work.The management of distributed energy sources and loads within the HRES network is facilitated through the operation of Local Controllers (LC) which are typically implemented as power electronic converters.Since the operation is controlled by electronic devices, the overall operation of the system is controlled by the Energy Dispatch Controller (EDC) or Micro-grid Controller (MGC) which sends command signals to every of the Local Controllers (LCs).The following subsections describe the mathematical modeling of each component of the proposed system to scrutinize its performance.

Solar PV
The apparoch proposed in this work takes into account two significant parameters, namely solar irradiance as well as ambient temperature, which could potentially affect the power output of solar PV at a certain period of time.The amount of electricity produced by the panels is calculated as a function of solar radiation 67 and the corresponding equation is given as follows 68 : where P pv is the output power (Watt) of the solar PV, P pv_r refers to the rated power (Watt) of a PV system under Standard Test Conditions (STC), G T denotes the hourly sun irradiance (W/m 2 ) that a solar PV panel surface receives.G T_STC refers to the solar irradiance at STC, which is equivalent to 1000 W/m 2 , ψ T denotes the coefficient of the temperature of solar PV (− 3.7 × 10 −3 °C−1 ), T amb symbolizes the ambient temperature (°C), T STC signifies the temperature of the PV cell at STC, which is typically 25 °C. (1)

Wind
Wind energy is a renewable and readily available resource that can be utilized to generate electricity.Electricity is generated by Wind Turbines (WTs) through the utilization of the kinetic energy of the wind 69 .The speed of the wind is the primary factor that governs the energy output of a WT.In order to obtain wind speed measurements at a specific hub height, it is necessary to convert the recorded wind speed at a different altitude.This is due to the fact that wind speed is not constant across different altitudes.The conversion process involves utilizing the power law equation of the wind profile, which is expressed in the following Eq.( 2) 70 .
The equation presented involves the variables V h (t) and V ref (t), which respectively represent the wind speed (measured in m/s) at the hub height H h (measured in m) of a wind turbine and the reference wind speed (measured in m/s) at the anemometer height H ref (measured in m).The symbol α indicates the friction coefficient, which is also known as the Hellmann exponent wind gradient, or power-law exponent.The variable α is subject to various parameters like as wind speed, temperature, the curvature of the terrain, and altitude above the earth's surface.The standard value of α is 1/7 for the areas with moderate roughness and good exposure 71 .The value of α can be calculated using the following Eq.( 3) 72 .
The output power P WTG (kW/m 2 ) from Wind Turbine Generators (WTGs) for a given wind speed V h (m/s) can be determined by using Eq. ( 4) 73 .
where P WTG_R is the rated power of WTG, and V h denotes the wind speed measured at the hub height of a WT.V R (m/s), V CI , and V CO symbolizes the rated wind speed, cut-in wind speed, and cut-out wind speed of the WT respectively.

Diesel generator
DGs have been identified as a dependable and cost-effective solution for mitigating energy storage requirements, particularly in rural industries as well as remote regions.DGs consist of compression ignition engines and an alternator that employs non-renewable resources specifically oil to generate electrical power.The utilization of the DG may vary depending on the system, as it can serve as either the primary or supplementary source.DG also serves as a backup source of energy in case of battery depletion during a period of high demand 74 .The quantity of fuel necessary for electricity production is dependent upon variables such as the generator's fuel heat www.nature.com/scientificreports/rate efficiency as well as the fuel's heat content.The efficiency of the generator is a function of the load present during operation.When developing a hybrid system, it is important to take into account the diesel generator's efficiency and fuel consumption per hour.The assessment of fuel consumption in a DG is based on the size of the DG (generator's power output) and the magnitude of the load it encounters.The majority of DGs function lies within the range of 80% to 100% of their designated power output.The fuel consumption (F DG (litre/kWh)) of a DG is computed by utilizing Eq. ( 5) 75 .
where P gen denotes operational or generated power output at time t and P R represents DG power rating of DG in kW respectively.The gradient and intercept coefficient of the fuel curve are both expressed in terms of litre/ kWh and denoted by α 1 and α 2 .The values of α 1 and α 2 have been determined to be 0.246 and 0.0842 respectively for the DG that is used in the design 76 .

Invertor
DC/AC and AC/DC power converters are essential for systems that utilize both alternating current (AC) and direct current (DC) elements.The solar PV panels and BES, deliver a DC output, however the load being considered demands an AC.The choice of the converter size is determined by the greatest demand for load (P PLD ).The equation for determining the inverter rating (P inv ) is as follows 77 : where η inv represents the efficiency of the inverter.

Battery energy storage system
Batteries serve the purpose of storing electrical energy for utilization during instances of power surges and the unavailability of alternative resources.Battery Energy Storage System (BESS) is commonly employed as a means of backup for hybrid stand-alone energy systems with the aim of enhancing their availability and facilitating load leveling for temporary fluctuations.Lead-acid batteries are frequently utilized as a viable alternative for HRESs due to their cost-effectiveness, low Depth of Discharge (DOD), elevated safety standards, recyclability, and rechargeable nature 78 .The capacity of the battery bank is dependent upon two factors such as the need for the electrical load and the projected duration of time at which the BESS is anticipated to provide power to the load.The BESS is expected to provide power to the load only when the RES output is inadequate to satisfy the energy demand.This duration of time is commonly referred to as the Days of Autonomy (DOA).The system's battery capacity C Bat in kW is determined based on the required demand and desired number of DOA, as specified by the following Eq.( 7) 79 .
where P L denotes the load demand (in kilowatts) that needs to be fulfilled by the BESS.D a represents the number of autonomy days, which is typically within the range of 3-5 days.DOD signifies the depth of discharge, which is set at 80%.Furthermore, η inv and η Bat corresponds to the inverter and battery efficiencies, which are 95% and 85%, respectively.
The process of charging as well as draining of a Battery Storage System (BSS) which primarily uses lead acid batteries occurs in a sequential manner.Consequently, the present state of charge of P BES (t) is reliant upon the prior state of charge of P BES (t−1).During the charging period, the updated state of BSS energy may be determined using Eq. ( 8) 80 .
where the symbol δ represents the self-discharge rate, η inv represents the inverter efficiency, and η BES represents the charge efficiency of the battery.When the output power of a RES is insufficient to meet the load demand, the BSS switches into a discharge mode.The updated status of BSS energy for the discharge phase can be determined by ( 9) 81 .

Case study
The approach implemented in this study can be used to create a compact independent PV-wind-battery-DG hybrid system which is shown in Fig. 3.This compact design is assumed to fulfill the energy needs of a small community located in a village and can serve as a potential case study location.The selected case study location for this research is Alleri village, situated in close proximity to Vellore city in the Tamil Nadu region of India.The geographical coordinates of the village are around 12.8095 latitude and 78.9913 longitudes.The solar irradiation, temperature, and wind speed data for the site location were collected from the database of NASA Power Data (5) Access viewer website database for the year 2021.The average value of solar irradiation and temperature during the said period are 5.17 (kW-h/m 2 /day) and 24.76 °C respectively.Also, the wind speed range is 5.06 m/s. Figure 4 depicts the hourly solar radiation, temperature, and Clearness Index data of the selected location and this is one of the data sets which is used to model the proposed HRES in this work.The main reason behind selecting the data set is that it has large variance over the period and it is useful to achieve realistic HRES design.This data is obtained from the NASA Power Data Access Viewer website 105 .The hourly wind speed statistic for the selected case study location is shown in Fig. 5.

Load profile
The load experienced by individuals is influenced by a variety of factors, including geographical conditions, economic circumstances, standard of living, as well as climatic conditions.In the context of usual household consumers, it is common for various electrical appliances to be utilized.These appliances include water pumps, electric irons, ceiling fans, LED lights, pedestal fans, televisions, refrigerators, etc.When assessing the energy consumption of these appliances, it is important to consider the average daily running hours for each appliance, as this provides a more accurate representation of their usage patterns within a household.The analysis of the load profile reveals that the load exhibits continuous variations over the course of the day.To enhance the authenticity or realistic nature of community load, a randomization factor of 10% can be added to the maximum value.The calculation of the total hourly as well as yearly load in kilowatt-hours (kWh) is dependent on the number of residential or domestic devices utilized along with their respective ratings.This calculation is performed using Eq. ( 10) 82 .
where E demand denotes the energy or electricity demand (kWh), n user represents the number of loads or equip- ment and T refers to the time of use.The examination of the load profile of a given area is of greatest significance when it comes to the development of a dependable and effective system customized to that specific region.The determination of battery size and modeling is dependent upon the load profile.Additionally, it is important to note that the reliability of the system is influenced by peak times and consumer behavior.These factors also play a significant role in determining the appropriate sizing of system components and the COE.The appropriate sizing of the HRES is mostly dependent upon the load requirements, and it signifies the energy consumption over one year.The data related to the electrical consumption of households of the selected case study on an hourly basis over a period of a year is obtained from the HOMER software.Figure 6 shows the daily load profile of the entire village.Figure 7 presents the annual hourly load curve.The daily rate of change in load is taken as 10%, whereas the hourly rate of change in load is considered as 20%.

Methodology Power management strategies
The integration of the Power Management Strategy (PMS) is a crucial aspect of hybrid energy systems.Its primary function is to ensure that a balance is maintained between the power generated by RES and the power required to meet the load demand.This is particularly important in situations where weather conditions are variable and there is a high degree of uncertainty associated with RES.As a result of the restricted amount of power that may be generated by renewable sources, the capacity of the DGs cannot rapidly be raised to satisfy the growing demand.In cases at which the power generated surpasses the demand, it becomes necessary to utilize a dump www.nature.com/scientificreports/load which is also known as a diversion load or dummy load to dissipate the surplus energy and protect the storage devices like battery banks from overcharging.The term "dump loads" is most commonly used to refer to deferrable loads, such as electric water heaters and hot air heaters used in residential environments, water pumps, and power resistances.Thus, it is essential to incorporate PMS in the design of such systems.In order to apply an EMS, it is necessary to conduct simulations that take into account the following conditions or approaches.State I: Battery charging approach: In case, at which overall energy output from RES surpasses the load demand within a specific time frame and the State of Charge (SOC) of the battery has not yet reached its maximum predefined capacity, the battery energy storage system is utilized to consume the excess energy for battery charging.Alternatively, if the aforementioned conditions are not met, approach 2 is initiated.
State II: Energy dumping approach: Analogous to the first approach, in cases at which the overall energy output of RESs exceeds the instantaneous load requirement and the SOC of the battery bank is at its maximum.Thus, in this scenario, the excess power is utilized in dump loads.The implementation of this particular strategy serves to mitigate the risk of battery bank damage resulting from overcharging.
State III: Battery discharging approach: RESs are unable to meet the load demand due to insufficient energy production.In this scenario, the foremost concern is to utilize the energy stored in the battery banks, instead of relying on the functioning of the DG.In this instance, the deficit in the production of electrical energy is compensated for using a battery.
State IV: DG supply approach: In the case that the energy produced by the combined RESs and the energy that has been accumulated in the battery fails to satisfy the energy requirements, the DG is activated in order to compensate for the energy deficit or to provide power to the load and to facilitate the recharging of the battery banks.The proposed dispatch strategies are represented in a simplified flowchart, as shown in Fig. 8.

Formulation of optimization problems
In order to develop a cost-effective and optimal performance of HRES, it is essential to prioritize the sizing of the various system components.The primary determinants are the ideal sizing of PV systems, WTGs, BESS, and

Objective function
The lifespan of a system can be significantly impacted by the integration of RESs and the utilization of components of superior quality.This in turn can lead to a reduction in the COE for end-users residing in remote regions.
Economic or cost assessment.Conducting an economic analysis is a significant consideration for any executing agency before the establishment of a renewable energy project.Prior to initiating a project, it is essential to conduct a comprehensive evaluation of its viability across various aspects including its economic feasibility.The economic feasibility of an HRES can be assessed through the evaluation of three parameters namely ASC, NPC, and COE.The Annualized System Cost (ASC) concept is used for economic evaluation.The solution with the lowest ASC is regarded as the most optimal provided it meets all other requirements and characteristics.The ASC can be described by the following Eq.( 11) 70 .
C PV is the cost of a solar PV panel per kilowatt, C WT is the cost of a wind turbine per kilowatt, C BES is the cost of a battery per unit, and C INV is the cost of an inverter per kilowatt.C INV represents the cost of the inverter, whereas P INV indicates the rating of the inverter.Each component's cost includes capital, replacement, operation and maintenance, fuel, and salvage expenses.
Furthermore, this research aims to decrease the total Net Present Cost (NPC) of the proposed HES in order to minimize the COE while maintaining optimal energy flow.The total NPC is a component that plays a significant role in the development of energy, and COE is the metric that is used for examining hybrid systems.The Total Net Present Cost (NPC) is a financial measure that can be employed to evaluate the economic feasibility of an investment project.The total net present cost (TNPC) can be determined based on the assessment of HRES 71 .The cumulative discounted cash flows throughout the year of the project's lifecycle are added together to determine the total NPC 79 .The total NPC may be determined using Eq. ( 12) 77 .
As pointed out earlier, the primary aim of this optimization problem is to minimize the COE by identifying the most cost-effective optimal size of the HRES component that can efficiently fulfill the energy demands.The present study incorporates the cost associated with the operation of HES, which involves the costs of PV array, WT, inverter, and DG.It is referred to as the cost per unit of electric power or the constant price per unit of energy.It is determined by employing the following Eq.( 13) 83 .
where the variable P load is utilized to indicate the hourly consumption of electricity or the aggregate load that has been served.The calculation of the NPC involves determining the present value of the total capital investment as well as operational costs incurred over the entire lifespan of the project.The parameter CRF is an abbreviation that stands for Capital Recovery Factor, which is a ratio employed to measure various components in relation to the current interest rate during a specific period which is given in Eq. ( 14) 33 .where the variables γ and τ represent the actual or real discount rate and the number of years.γ ′ represents the nominal interest rate and f denotes the annual inflation rate respectively.
Reliability assessment.The assessment of system reliability is of greatest significance as unpredictable weather can significantly affect the performance of power generation equipment.An effective power supply system necessitates minimal or ideally zero loss of load.The Loss of Power Supply Probability (LPSP) is a statistical metric that determines the possibility of power supply failure.This failure can be attributed to either insufficient renewable resources or technical issues that prevent meeting the demand for power.LPSP serves as an indicator of the reliability of the hybrid systems in fulfilling the power demand.It can be defined as follows 84 .
where t f represents the final time or operating period which is typically 8760 h for one year, P Load indicates total or overall electric energy demand and DPS is the deficiency of power supply.Chronological simulation and probabilistic approaches are the two ways to compute LPSP.The first method makes utilization of time-series data over a certain period, while the second relies on the energy accumulative impact of an energy storage device.LPSP has values ranging from 0 to 1.When the power generated equals the load demand, the load is fully supplied by ( 14) the generated power and the LPSP value is 0. The LPSP value of 1 indicates that the load is fully unsatisfied i.e. load demand is not met by the generated power.
System constraints.This study introduces the concept of the Renewable Factor (RF) as a means of evaluating and comparing the energy output derived from renewable sources with respect to the energy produced by DGs.Similarly PV fraction, WT fraction also measures the respective energy output aganist energy produced by DGs.The ideal system is one that only uses renewable resources which is demonstrated by the variable RF and in this case the value is 100%.Although, the RF = 0% demonstrates that the energy produced by a DG is equal to the energy produced by renewable resources.RF, PV fraction (PVF) as well as WT fraction (WTF) are determined by employing Eqs. ( 17), ( 18) and ( 19) 85 .
where,N max PV ,N min PV , N min WTG , N max WTG are the number of maximum and minimum of PV and WTG respectively.
Social parameter assessment.The Human Development Index (HDI y,z ) is commonly utilized as a statistical tool for assessing a country's overall growth in its economic as well as social indicators.The HDI is determined by combining the values of the elements "y" along with the control approach "z" using a specific Eq.( 22) 86 .The HDI is a numerical metric used to assess the level of continuous advancement in a country or region.It relates to the overall well-being of persons which includes the services, physical infrastructure, and resources available to them.This research incorporates the HDI into the objective function to provide a design solution that considers the well-being of society.Excess or wasted energy may improve the overall quality of life in society.The HDI is calculated in the following manner 87 .
where F max_DP represents the percentage of the greatest energy that is dissipated.F max_load is the limiting factor for the extra load that the available power can handle, while ℇ pop indicates the number of persons benefiting from the hybrid system.

Proposed optimization technique
The optimization of hybrid systems through the implementation of an EMS can be classified into four main categories.These categories include load supply optimization, which is a simple optimization method dependent on load as well as power.Technical optimization is another category that aims to improve the operation and durability of the hybrid system.Economic optimization on the other hand focuses on reducing maintenance costs through an optimization algorithm based on cost reduction function.Lastly, techno-economic optimization is a more complex category that employs an optimization algorithm depending on the objective function.As already said, this research employs a 3E analysis to evaluate the energy, ecological, and economic aspects with social, and technical factors assessment of HRES.To address these issues, intelligent methods like MH algorithms are the best way to find a solution.Figure 9 illustrates the problem-solving methodology adopted in this work using the mHHO optimization technique.The flow chart (Fig. 9) depicts the usage of real-time input datasets to evaluate objective functions, namely technical, economic, environmental, and social factors.The primary purpose is to achieve output segments that include the least COE, NPC, and ASC, minimum LPSP, and maximum reliability of the system while maximizing the HDI.

Harris Hawk optimization algorithm (HHO)
The HHO algorithm takes inspiration from the hunting as well as capturing strategies used by hawks.The objective of this method is to initiate an assault on a target, often a rabbit, by the coordinated efforts of many hawks approaching from various angles, with the intention of catching the prey unaware and unprepared.The typical HHO method may be subdivided into two fundamental stages, namely exploration and exploitation.The exploitation stage consists of two distinct stages, namely the hard passage and the soft passage 50  www.nature.com/scientificreports/Exploration Exploration is often regarded as a global search mechanism within the context of optimization algorithms.During this stage, the algorithm systematically investigates the search space in order to identify and evaluate potential solutions that exhibit desirable qualities.In the context of HHO, the Harris hawks exhibit a behavior of perching in different places without a specific pattern, using two distinct approaches to identify potential prey or rabbit.Consider that q is a probability for both the exploration in randomly placed tall trees as well as the exploration based on the location of the rabbit and other hawks as indicated in the following Eq.( 23) 88 .
where t represents the number of iterations.X(t) denotes the current location of the Hawk, while X(t + 1) rep- resenting the position of the Hawk in the subsequent iteration.Additionally,X prey signifies the current position of the rabbit.The variable X rand (t) represents a stochastic location for the hawk.The variable X m (t) represents the average location of the hawk's current location, which is determined using Eq. ( 24).The lower and upper limits of variables are denoted as LB and UB.The factors r 1 ,r 2 ,r 3 ,r 4 and q are assigned random values.The haw's position is given by the following equation where X i (t) specifies each hawk's position in iteration t, and N stands for the overall number of hawks.The transition from exploration to exploitation is dependent upon the escaping energy of prey.In order to represent the process of transitioning from exploration to exploitation, the energy level of prey is conceptualized and represented in a mathematical model as follows.
where E represents the prey's escaping energy, E 1 denotes the linearly decreasing factor with a range between 2 to 0, T denotes the maximum number of iterations, and E 0 signifies the beginning state of the prey's energy.Exploration occurs when |E| ≥ 1 but exploitation takes place in subsequent stages when |E| < 1.

Exploitation
Exploitation is often seen as a localized search in comparison to the exploration stage.The algorithm attempts to minimize big jumps inside the search space while simultaneously refining the outcome during this particular phase.There are four internal phases to the exploitation phase, which are based on the level of escaping energy (E) and the probability (r) associated with the successful escape of the prey (r).

• Soft besiege
During this particular stage, when |E| ≥ 0.5 and r ≥ 0.5 the prey has sufficient energy and attempts to escape by unpredictable and deceptive movements.However, ultimately it is unable to successfully escape.During these efforts, the Harris' hawks perform a strategic movement of circling the rabbit softly in a gradual ( 23) www.nature.com/scientificreports/manner, aiming at causing exhaustion before executing an unexpected attack.This behavior is represented by using the following Eq.( 27).
where �X(t) represents the difference among the position vector of the prey and its current position at iteration t.Additionally, r 5 denotes a randomly generated value that falls within the interval (0,1).The variable J denotes the random jumps made by the prey all over the process of escape.

• Hard besiege
During this particular stage, the prey's energy is reduced to a point where it is unable to escape, denoted by |E| < 0.5, r > 0.5 .Consequently, the hawk employs a strategy including hard circles in order to capture the prey, as shown in Eq. ( 30).

• Soft besiege with progressive rapid dives
The prey has sufficient energy to effectively escape when |E| ≥ 0.5 but r < 0.5 , and a mild besiege is still built before the surprise attack.In comparison to the prior situation, this technique is more sensible.Observing hawk behaviors, it is hypothesized that they may choose the optimal dive to grab prey in competitive scenarios.Thus, to softly besiege, it is assumed that the hawks could assess (decide) their next step using the following Eq.( 31) However, as a result of the deceptive and unpredictable movement pattern shown by the prey, the Hawk adjusts its present location to position (Z) by following Eq.(32).
where D denotes the dimension of the problem, the random vector is represented as S with dimensions (1 × D) , and the levy flight function is denoted as LF as stated by the following Eq.(33).
where u and v represent random variables that are uniformly distributed between 0 and 1.The constant β is assigned a default value of 1.5.Therefore, the ultimate approach of updating the locations of hawks during the soft besiege phase may be executed using the following Eq.( 35).

• Hard besiege with progressive rapid dives
During this phase, the prey may have a limited chance of escape due to a significant fall in its energy, denoted as |E| < 0.5 , and r < 0.5.Consequently, an effective besiege is established prior to executing a sudden attack in order to capture as well as kill the prey.The hawk reduces the circle of its hunting area and adjusts its locations by (35).The locations denoted as (Z) and (Y) are defined by using Eqs.( 32) and (36), respectively.With respect to the described methodology above, there exist several factors that require optimization in order to enhance performance.In the modified method which is described in the next section, the escaping energy of the prey is optimized with 9 different operators which are called as inertial weights.The detailed description of energy change with respect to nine different operators is described below.( 27) The present work introduces a new approach, namely the Modified Harris Hawk Optimization Algorithm (mHHO), to address the sizing as well as the techno-economic assessment of HRES.The mHHO algorithm is capable of transitioning from a phase of exploration to an exploitation phase, and thereafter adapting its exploitative behaviors dependent on the energy level of the prey's attempts to escape or the prey's escaping energy.In the context of the HHO algorithm, the concept of energy is often linked to the fitness or objective function which undergoes optimization.The search process may experience a decreasing exploration-exploitation balance as the energy component decreases.According to Gauri Sahoo et al. 89 , the transition from the exploration to the exploitation phase is in the context of the prey's escape energy value inside the HHO algorithm.Excessive exploration of the search region might bring a degree of unpredictability or randomness and could lead to being stuck in local minima.In addition, excessive exploitation may lead to reduced levels of unpredictability and could prevent the achievement of optimal outcomes.Hence, an appropriate balance needs to be conserved between the exploration and exploitation phases during the search process 50 .In order to achieve the balance in the current work, various operators are being primarily used for updating the position.In order to symbolize the decreasing energy factor in the HHO system and to optimize the escaping energy one can make use of these operators to generate mHHO algorithms.In this work, the HHO algorithm is modified using 9 different operators as shown in Table 2.The mathematical model or the formula of the operators is provided to further know its behaviors in optimizing the escaping energy.Incorporating the formulae given for various operators in Table 2 on Eq. ( 24), the new escaping energy of the prey can be evaluated.Once the expected energy change is obtained using the respective operators, the performance of the algorithm is evaluated effectively.
When incorporating these different operators on HHO, the transition from exploration to exploitation phase in the HHO algorithm is achieved.This further helps to achieve optimal solutions.The optimal solution and functioning of the modified HHO (mHHO) using these 9 different operators are tested on 10 different problems proposed by CEC 2019 to demonstrate or show its effectiveness.Moreover, from the implementation of the 9 operators on CEC 2019, one can select the best operators for sizing and techno-economic analysis of HRES.

Results and discussion
This section presents a detailed examination of the efficacy of the proposed mHHO algorithm and the performance assessment of the HRES system employing mHHO and other MH algorithms.This section is divided into two parts.The first subsection focuses on the proposed Modified Harris Hawks Optimization (mHHO) algorithms.It performs a comparative study among nine different versions of the algorithm using CEC 2019 benchmark functions.This analysis helps identify the optimal mHHO version.Subsequently, this optimal mHHO is compared against various other MH algorithms.The second subsection then evaluates the effectiveness and robustness of the chosen mHHO algorithm.Here, the algorithm is applied to a Hybrid Renewable Energy System (HRES) optimization problem and its performance is compared with other MH algorithms such as Chinese Pangolin Optimizer (CPO) 99 , Arithmetic Optimization Algorithm (AOA) 100 , Exponential distribution optimizer (EDO) 101 , Liver Cancer Algorithm (LCA) 102 , Novel Bat Algorithm (NBA) 103 , and HHO 50 .The parameter settings for every approach are obtained from the standard sources and are shown in Table 3 for ease of analysis and reference.All the simulations are performed by using MATLAB: 2021a on VivoBook_Asus having an Intel core i5 processor (Intel(R) Core(TM) i5-1035G1), 8 GB RAM, and Windows 11 ultimate operating system.

mHHO performance analysis
In order to evaluate the effectiveness of the proposed algorithm, the mHHO has been implemented on a total of 9 distinct benchmark test functions 104 .To substantiate the performance, the statistical tests namely Friedman's rank test (f-rank test) and convergence profiles, are used to demonstrate the improved performance of the proposed algorithm.It is important to note that the calculations have been derived from a total of 51 runs, ensuring a robust and reliable assessment of the algorithms' characteristics.For fare comparison, maximum iteration and population size of all these algorithms are set to 500 and 50 respectively.The f-rank test includes Table 2. Different operators used for improving the HHO algorithm.

S. no.
Operators or factor Formula References 1 Linear decreasing factor assigning a distinct rank to each algorithm being compared based on its performance.The results of the f-rank for each benchmark function are presented in Table 4 for each function separately.The overall total f-rank has been determined by considering the performance across all test functions and is also shown in Table 4. Also, all mHHO are compared in terms of mean and standard deviation values.After observing the results in Table 4, the performance of HHO_NatExpo is superior on all benchmark functions, since it achieved the highest overall ranking among all the tested algorithms.According to the results presented in Table 4, HHO-LogDec produces superior performance compared to other cases when applied to functions F1, F4, as well as F9.For function F3, HHO-Oscilla inertia weight shows better results as compared with the other algorithms.HHO-SigDec shows superior results as compared to the other operators' results for the function F5.For functions F6 and F8, the HHO-Chaotic operator performs better as compared to others.In the context of function F7, it has been observed that the utilization of HHO-LinDec produces superior performance in comparison to other approaches.HHO-ChaoRand outperforms competitors for functions F10.The results are also validated for each case using an f-rank statistical test.
From the overall f-rank value and average f-rank, it is evident that HHO-NatExpo is showing remarkable superiority over the other approaches and achieves the 1st rank followed by HHO-ChaoRand.In addition to the obtained simulation results, the convergence profiles of the mHHO algorithm using all operators are shown for a total of 500 iterations.The convergence characteristics for all benchmark functions (F1-F10) are also shown in Fig. 10.Interestingly it has been observed that the convergence profile of the HHO-ChaoRand and HHO-NatExpo approach is superior on all benchmark functions compared to other modified algorithms.Additionally, Fig. 11 displays a box plot analysis to further prove the performance of the algorithms on the benchmark functions.Based on the analysis of box plots, it is clear that the mHHO algorithm utilizing the Natural Exponent Operator (HHO-NatExpo) exhibits a lower median value and also a smaller interquartile range.
Also, from Table 4 it is clearly observed that the HHO_NatExpo secured the First rank.This result necessitates the selection of the HHO_NatExpo algorithm for comparison with the other MH algorithms for further designing of HRES.The outcomes presented in Table 5 show that the HHO_NatExpo can outperform other algorithms for functions F1, F2, F3, F4, F8, F9, and F10.In the sixth function, the EWOA algorithm gives the best possible outcomes in terms of mean and standard deviation values.For functions F5 and F7, the PSO algorithm shows a superior result as compared to the other algorithms.As a result, the HHO_NatExpo algorithm is shown to be successful for seven out of a total of ten numerical test functions, while the EWOA algorithm is effective for one function, and the PSO algorithm is effective for three functions.Overall, it has been determined that the HHO_NatExpo algorithm is the most successful algorithm for the CEC 2019 benchmark function.

Statistical testing
Two non-parametric tests Wilcoxon's rank-sum (p-rank) and Friedman rank (f-rank) test are performed to confirm the superior performance of the HHO_NatExpo algorithm.For analyzing the outcomes, wins (w), losses (l), and ties (t) are the first step in calculating the p-rank.The algorithm that wins against the HHO_NatExpo is assigned with a " + " sign to show its superior performance.Similarly, if the algorithms show worse outcomes as compared to the HHO_NatExpo algorithm, then assigned with the "−" sign.Finally " = " sign is assigned when there is no statistical difference between the compared algorithms.The HHO_NatExpo algorithm outperforms other MH algorithms in most of the functions and one can observe easily the result in the W/L/T row in Table 5.Additionally, the f-rank test assigns a distinct f-rank to every algorithm being evaluated and it is tabulated in the fourth row of every function in Table 5.The average f-rank is calculated for each algorithm and it is found that the proposed HHO_NatExpo algorithm achieved the Ist rank followed by the PSO and SCA algorithms.

Results and discussion of HRES
The effectiveness and robustness of the proposed mHHO (HHO_NatExpo) algorithm on HRES is illustrated in this section since it is found to be optimal and the best one for designing reliable systems in comparison with other mHHO and MH algorithms.The simulation is performed with the population size set to 30 and the maximum iteration considered as 250.Furthermore, the technical and economic specifications for the simulation of HRES for components like Solar PV, Battery, DG, and others are presented in Table 6.The project's duration is taken as 20 years, while the interest rate is fixed at 3.5%.Table 3.Details of various algorithm parameters.The MATLAB 2021a program is used to simulate the experimental data.The simulation time step and duration considered are 1 h and 1 year respectively.The maximum ratings for solar PV panels, WTs, and number of batteries have been set at 300, 300, and 170 respectively for the purpose of comparing results.The optimum results includes the maximum rating of solar PV panels, WTs, the total quantity of batteries, and the highest possible rating for the DG.The viable and optimum solution is evaluated based on ASC, NPC, and COE.Table 7 displays the entire set of optimized results achieved for the case study using the proposed mHHO along with other MH algorithms such as basic HHO, LCA, MPO, AOA, EDO, NBA, and SCA algorithms.The findings proved that the mHHO algorithm yielded the lowest NPC and ASC, resulting in the least COE.The mHHO algorithm estimates the setup of a 299 kW solar PV system, 299 kW wind turbines, 169 batteries, and a 94 kW DG with an Annual System Cost (ASC) of $ 1,16,090 resulting in a COE of 0.1413$/kWh.From Table 7, it is observed that mHHO produces an 89% decrease in the value of COE against LCA, and also with basic HHO it offers 0.07%.Also, in comparison with other MH algorithms such as MPO, AOA, EDO, NBA, SCA, the proposed mHHO achieves around 0.49% decrease in the value of COE.This proved to be one of the most economical designs considering the ground realities present at the site location.The COE obtained using mHHO indicates that the proposed system can provide electricity to off-grid locations at a reasonable cost.Furthermore, Fig. 12 displays the monthly fluctuations in load demand along with the different components of power production of the HRES system over the course of a year.It can be seen from Fig. 12 that RES generates significant power during June and July making the period as one of the high-power production months.On the other hand, October is the least productive month in terms of power generation leading to increased utilization of battery power and DG.

Algorithm Parameters
Table 8 displays a comprehensive annualized cost evaluation of the hybrid system.The annualized costs of the system are determined using Capital Recovery Factor (CRF).Table 8 provides a thorough analysis of the cost of individual components and the annualized cost for the most efficient hybrid system.Yearly expenses are calculated using the CRF. Figure 13 demonstrates a comparative comparison of the cost of different components of HRES using various techniques.The results in Table 8 and Fig. 13 illustrates a significant decrease in the cost of several components when using mHHO leading to a reduction in the total annualized cost.Figure 14 shows an overview of the convergence characteristics of all algorithms.From an enhanced view of Fig. 14, it is observed that both the HHO and mHHO algorithm performs better compared to other MH algorithms but after the 165th iteration mHHO algorithm is giving better results than the basic HHO.Compared to the mHHO method, other algorithms required a longer time to simulate the proposed system.The ASC, cost of individual components, convergence curves, and sizing make the mHHO an imperative optimal algorithm for further HRES design.
It is observed that the proposed approach meets the whole energy requirement using solely solar, wind, and batteries except a infinitesimal duration using DGs. Figure 15 illustrates the monthly average energy balance and usage during a period of one year.It should be highlighted that the amount of energy produced by solar and wind has matched or is equal to the available resources (solar and wind).In October, when wind turbines are producing less power, a DG is used to meet the electricity need.Moreover it is revealed that during a given month, DG will be in operation whenever the cummulative energy produced by solar, wind and battery fails to satisfy the load demand.Throughout the remaining months in a year, a significant quantity of wind energy is produced because of the greater accessibility of natural resources.Except summer season, the utilization of www.nature.com/scientificreports/battery banks is increased which means the power extracted from batteries is greater.It can be also seen from Fig. 15 that only in 3 months (June, July, and August), where more excess energy is available compared to other months of the year.Notably, wind energy makes up a significant portion of the contribution.Table 9 provides a concise comparison of the energy generated by every element of the designs derived from HHO, mHHO, LCA, EDO, MPO, NBA, and SCA.The total extra energy for the whole year utilizing mHHO is 209,100 kWh/year as shown in Table 9, which is lower than the results obtained using other techniques.
To better comprehend the power exchange among all the different elements of the system in a better way, a total of 2 weeks have been chosen.One week is taken in February when the load is lower whereas the second week is in October when the demand is greater.DG is used at these hours to supply electrical power to the system.Figure 16 displays a power exchange during 1 week in October.During this week DG is used when there is insufficient power from solar, and wind sources and also the battery's SOC is at or below the minimum threshold of 20% as outlined in the PMS.One can easily observe from Fig. 16 that, out of 168 h, the DG work for only 7 h in a week using mHHO.This has translated to 96% working of renewable sources with batteries by reducing the CO 2 emissions to a bare minimum value.
Similarly, Fig. 17 illustrates the energy management during the last week of February.The data in Fig. 17 indicates that sufficient solar power and wind power are produced as a result of the abundant availability of both renewable sources.DG power is not required since the complete load requirement is satisfied only by batteries, solar, as well as wind energy sources.In addition to energy analysis, the total energy production and load characteristics for the complete year are displayed in Fig. 18.From Fig. 18, it is observed that the total output power (kW) surpasses the load (kW) showing that the total output power is enough to meet the load requirements.
Battery SOC evaluation is a crucial concern in systems that use batteries for storage.Figure 19a,b displays the fluctuations in the battery bank's SOC along with the input and output energy over the course of a year.The starting SOC and the lowest acceptable SOC are set at 100% and 20% respectively.Figure 19a shows that the battery SOC consistently stays within a predetermined limit.To understand better, it is considered that on January 1st, at 0000 h, the initial SOC of the battery banks is 100%.One can observe clearly from Fig. 19a that the SOC dropped to a minimum of 20% for only a short time during the whole year.Figure 19a demonstrates that the battery SOC is generally optimal with some exceptions during periods of limited natural resources.For example, during the month of January and high load demand in October-November, the SOC drops from its maximum value to a certain limit.This corroborates the result of minimal DG operation which is reported also in Fig. 16.If there is excess energy available above the charging rate and if SOC is also 100%, the excess energy may either be stored as a delayed load or dumped load.If the energy discharge exceeds the battery discharge rates, a DG will be used as the power generation source.
Table 10 presents a comparative assessment of technical, social factors, and environmental factors such as LPSP, HDI, RF, WTF, PVF, and CO 2 emissions using different MH algorithms for the proposed design.Based on the results, it is observed that mHHO outperforms other MH algorithms.Further, the comparative pictorial representation of these factors are displayed in Figs.20 and 21.The CO 2 emissions of the proposed HRES system have been reduced from 6165 kg/year using LCA to 322.7402 kg/year as reported by the mHHO method.The CO 2 reduction is 94% which is a remarkable reduction.Also, in comparison with the basic HHO algorithm, the mHHO shows better results with a 0.9476% reduction in CO 2 emission.The proposed mHHO also shows prominent results in HDI.It produces the HDI of 0.6550 which is almost advocated by WHO standards.From the above analysis and results, it is clearly proved that the mHHO is one the strong algorithm for implementation of HRES.

Conclusion
This work presented a new approach known as mHHO optimization algorithm for the 3E analysis of HRES.Primarily the paper is focussed on Energy assessment, Economic evolution (COE, NPC, and ASC), and Environmental analysis (PVF, WTF, RF, CO 2 emission).In addition, the technical and social factors such as LPSP and HDI for HRES (PV-Wind-Battery-DG) are also evaluated.To obtain a practical HRES, mHHO optimization algorithm is utilized to achieve the optimal system design and to determine the appropriate sizes of HRES components.The proposed mHHO algorithm is modified using 9 different operators and tested on CEC 2019 benchmark test functions.The comparative analysis is carried out based on statistical tests such as the f-rank test with convergence and box plot studies.Then, the algorithm that has yielded a higher rank i.e. mHHO (HHO-NatExpo) which has good performance and a good convergence profile compared to others is selected for the HRES problem solving in the chosen case study location.The proposed mHHO algorithm demonstrated superior performance by achieving the lowest COE of 0.14130$/kWh, ASC of 116,090$/year, and NPC of 1,649,900$.Furthermore, mHHO has exceptional efficacy in reducing the environmental pollutant CO 2 emission of which is accounted to a mere 322.7402 kg/year while maintaining HDI at the prescribed level.From an optimization perspective, it has been demonstrated that the mHHO algorithm has superior outcomes in 3E assessments compared to basic HHO and other MH-algorithms (LCA, MPO, AOA, EDO, NBA, and SCA algorithm).The results of the simulation indicated that using renewable resource might potentially serve as high reliability for optimizing hybrid energy systems.The incorporation of RES into the hybrid system has been seen to provide the advantage of meeting the entire energy demand by reaching 100% fulfillment in the near future.On the other hand, one of the crucial problems to address is the examination of the mHHO algorithm's behavior throughout the exploration and exploitation phases.Due to the algorithm's stochastic nature, there is a significant probability that it can diverge from the global optimal point and get stuck in local optima thereby reducing its reliability for different optimization challenges.In addition to said limitations, further the incorporation of adaptive parameters does not consistently provide a well-balanced functioning and may sometimes lead to premature convergence.
Based on the findings and implementation challenges, the following recommendations are suggested for future work enhancement and development.A laboratory-scale experimental setup can be developed to facilitate deeper research into the performance of HRES in harsh environments.The simulation model may also be created to understand the HES's power flow regulation process.Consequently, this work helps to raise potential awareness about the installation of RESs in different parts of the rural areas intending to facilitate an economically sustainable energy supply. https://doi.org/10.1038/s41598-024-70663-5

( 11 )
ASC = f (N PV C PV + N WT C WT + N BES C BES + P INV C INV + N DG C DG )

Fig. 7 .
Fig. 7. Yearly load profile of the entire village.

Fig. 12 .
Fig. 12.Output Power of all components and load demand.

Fig. 13 .
Fig. 13.Cost of different components of HRES system using various algorithm.

Fig. 15 .
Fig. 15.Assessment of monthly energy usage for the proposed case study.

Fig. 16 .
Fig. 16.Energy balance with SOC for the 3rd week of October.

Fig. 17 .
Fig. 17.Energy balance and SOC of battery for the 1st week of February.

Fig. 19 .
Fig. 19.(a) State of charge of battery (%).(b) Total energy input and output of the battery during the period of 1 year. .

Table 4 .
Comparison of simulation outcomes of HHO algorithm with different inertia weights (CEC 2019 benchmark functions).

Table 5 .
Comparison of simulation outcomes of HHO algorithm with different algorithms (CEC 2019 benchmark functions).

Table 6 .
Technical specifications and associated expenses for various components of the HRES.

Table 7 .
The optimal sizing results obtained using different algorithms.

Table 8 .
Cost of each component of the proposed system by using different algorithms (PV-WT-BES-DG).

Table 9 .
Assessment of energy for various algorithms for the proposed system (PV-WT-BES-DG).

Table 10 .
Comparison of technical and social factors for the proposed system using several algorithms.